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Abstract 



Using the polarizability of a free electron gas in a magnetic field and the 
Current-Density Functional Theory (CDFT) developed by Vignale and Ra- 
solt, we derive the gradient and current corrections for the energy functional 
of a non-uniform electronic system in a strong magnetic field. 

First, we find the Tomishima-Shinjo functional by neglecting the current 
variation. Taking into account the current variation leads to new gradient 
terms which change contraction effects in the direction perpendicular to the 
magnetic field. 
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I. INTRODUCTION 



In a magnetic field the electronic density is elongated along the field direction. There 
are many intuitive ways of seeing this. For example, in a path integral description, since 
the zero field propagator is multiplied by the phase factor exp[j^ § A ■ dl], each closed path 
is reweighted by a phase proportional to the magnetic flux it encloses. In the average over 
paths which produces the electronic density, paths parallel to the field, enclosing no flux, 
receive a weight of one; while paths perpendicular to the field, enclosing flux, receive a 
reduced weight. The resulting average for the electronic density, thus, has an elongated 
shape. 

The lowest order statistical model however leads to an isotropic densityfl]. The 
anisotropy requires gradient corrections [2]. The necessity of going beyond the lowest order 
has also spawned a variety of approximate methods to regain anisotropy [3,4,5,6,7,8]. 

In this Letter, we first review the fundamental results of CDFT. Then, neglecting the 
current variation, we rederive the gradient corrections of Tomishima and Shinjo. Finally, 
taking this variation into account, we calculate all second order electronic density gradient 
corrections. 

II. CURRENT DENSITY FUNCTIONAL FORMALISM 

For systems in a magnetic field, Vignale and Rasolt [9] show that the density p and the 
paramagnetic current density j p , related to the gauge-invariant total current density, by 



uniquely determine the external scalar, V x (r), and vector potential A(r). The ground state 
energy functional, which may be written as 



j(r) =j p ( r ) + — P (t)A(t) 



(1) 




(2) 



2 



where K[p,j p ] is the kinetic energy of a free electron gas in the external potentials V x (r) 
and A(r), and E xc \p,j p ] is the exchange-correlation energy, is minimized at the ground 
state energy by the correct density and current distribution. In the spirit of Thomas-Fermi 
approximations we subsequently ignore the exchange-correlation term and the Hartree mean 
field term. 

Although K[p, j p ] is not invariant under a gauge transformation: 

A (r) - A(r) - V X 

j p (r)^j p (r) + ^p(r)Vx (3) 

the difference: 

UWh\ = KWh\ -jj dr ^y = E ^h\ - J drp(r)V x (v) " f (4) 

which depends on the energy and physical current, is gauge invariant and, since unchanged 
by the addition of the gradient of an arbitrary function, depends on j p only through 
u{r) = V x (jp/p), the vorticity introduced by Vignale and Rasolt [9]. 
In the next two sections we determine the gradient expansion for U . 

III. GRADIENT CORRECTIONS IN STRONG FIELDS 

We consider an electron gas in an uniform magnetic field B = B z. As a first step we 
will neglect the variation of vorticity. Then the functional U is reduced to a functional of 
only one variable, p. We postulate an expression for U with gradient terms: 

U[p] = J rfr{/(p(r))+^(p(r))|V ± p| 2 + ^||(p(r))|V|| P | 2 } (5) 

where _L and || denote orientations perpendicular and parallel to the field. 
The function f(p(r)) can be determined from the Euler equation 

= p . (6) 

Sp(r) 



For the assumed U this gives 



f(p(r))-2^( /3 (r))Vip-2^||(p(r))Vjp 

- 2l(p(r))|Vxp| 2 - sf|(p(r))|V|| P | 2 = -K(r) + p 

which, at uniform density po and no external potential, is simplified to: 

f'lfto] = Po 



(7) 



(8) 



where po is the chemical potential of a free electron gas. 

<?_i_(p(r)) and g\dp(r)) will be determined by comparing the calculated linear density 
response to a variation in the external potential, <5p(q) = x( ( l)^^x(q)9 to the free electron 
gas results. A second variation of the Euler equation gives (after a Fourier transform), 



/"(A)) + 2g ± (p )kl + 2g l{ (p )k 2 = -l/x(k) 



(9) 



at uniform density po- 

As a preliminary example recall the case of no magnetic field. The chemical potential 
p = h 2 kp/2m = h 2 (3% 2 p) 2 ^ 3 /2m so eqn. (§) reproduces the usual Thomas- Fermi result, 



f(p) = ^-(3vr 2 ) 2 V/ 3 
10 m 



(10) 



and the Lindhard response formula 



X(k) 



mkp 



1 2(A; 2 - k 2 



In 



2k f 



2k F + k 



ink i 



h 2 n 2 V 3' 



\-hk/2k F f) +0(k 4 ) (11) 



or 



X-\k) 



2^2 



mkf 



l + l -(k/2k F ) 2 ^+0{k*). 



(12) 



gives, on comparison with eqn. @, the von Weizsacker gradient corrections 

^ ± ^" 3m 8A;|, 9 8pm 
(with the 1/9 factor from the long wavelength comparison) [10]. 



(13) 



In the case of very strong fields, where all the electrons are in the lowest Landau level 
and spin polarized, the free electron gas chemical potential p = hu c /2 + h 2 k 2 F /2m with 
k F = 2n 2 Tip/muj c so 

The free electron gas susceptibility for this case is well known since Horing [11,12] 

2k F + k» 



X(k) = -AnWh k\\ 111 



(15) 



2k F — k\\ 

where F(k) = e ~ fc2 ' 2 / 2 and / is the associated magnetic length defined by: I 2 = hc/eB . The 
expansion of x _1 to second order in k gives: 

,„ , Arr^p 2tt 4 ^ 5 p 1 2 ^ , 2 

-* (k) = ^7 ki " 12^" (16) 

so 

^tfp 

9± = — 4-3 17 ) 

and 

9\\ = - — — (18) 

The resulting energy functional: 
is the same as the one originally derived by Tomishima and Shinjo [6]. 



Note that even though the free electron states are still plane waves along the field direc- 
tion the gradient correction for this direction has the opposite sign from the free field case. 
This is a consequence of the 1-dimensional nature of the high field limit. Explicitly, the free 
fermion response function in 1-dimension is 
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X(k) 



-m 



h z irk 



In 



2k F + k 



2k F - k 



or 

;2 



X -\k) = {l - ^A; 2 /^} + 0(A; 4 ) (21) 

with kp = vrp. Therefore [13] 

7T ft 2 1 fi 2 

^" ~ ~24J^m ~ ~24p~m ^ ^ 

which is exactly the high magnetic field result derived above. 

IV. CURRENT CORRECTIONS 

The Tomishima and Shinjo approximation is obtained when the variation of vorticity is 
omitted. However, this variation involves the gradient of p and can not a priori be neglected. 
We therefore generalize our previous expression for U to include vorticity as follows: 

U[pM = J rfr{/(p(r),z/(r))+^(p(r),Kr))|V ± p| 2 + ^||(p(r),z/(r))|V,| P | 2 } (23) 

From the definition of the vorticity: 

i/(r) = A + Vx (±). (24) 
mc p 

Consequently, we do not need to add in the U[p, u] expansion terms in | W| which are higher 
order. 

Variation of the energy with respect to j p , 5E/5j p = 0, using the chain rule result 

5 15 
-V x 



<5j P (r) P( r ) M r ) 
gives: 



j = -- Vx l^ ( 25 ) 
m ov 



to second order. The Euler equation, 5E/5p(r) = p, becomes: 
df 

^[p, v]-2 gi ylp-2g\\V\p 



^|V,p| 2 -f|V„p|^-K( r ) + p + ^ (26) 
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At uniform density we get: 

?L 

dp 



[Po, ^o] = p[Po, M with u = u c z 



so the expression for / is formally the same as ( |i~4|) with the variables p and v: 

h, , 2vr 4 ^ 4 a 



f[p, 



3m 3 1 



IP 



(27) 



With equation 



we can now calculate the total current density, 

.4*4 9 



ft 47T 4 ^ 4 

1 2m 



■J>Vp x u 2 



(28) 



m l^o I 

where we have replaced \v\ by |^o| in the second member of (|28|) since the z/ variations 
are higher order. In the limit of very strong fields we find the result of Skudlarski and 
Vignale [14]. On the other hand equation @ remains unchanged at po since all the additional 
terms are proportional to Vp. Thus, the expressions for g±, g\\ are absolutely the same as 
in section III with uj c — ► \v\. 

In return, E is modified by the presence of the current, and for a comparison with 



Tomishima and Shinjo we have to develop u, eqn. pi), using eqn. fl28|): 



h 

2m 



N±P\ Vip 



4tt 4 ^ 4 
m 4 ujr 3 



|v ± p| 2 + pvlp 



(29) 



Using this expansion in eqn. (|27|) for f[p, u] and with eqn. (j28[) for the current we finally 
obtain for the energy functional E: 



E[p}=Idr[^p + 



lr |V ± p| 2 3~'// 5 ^ , 



3m 3 LU r 



im p 



+ 



m 4 u; c 3 



5m 



' 7T 4 ^ 4 ' 

m 4 a; r 3 



P 3 |V ± p| s 



ft 2 \V\\ P \ 
24m p 



dr p!4 



(30) 



V. CONCLUSION 



Current corrections do not change anything along the field direction since this correction 
is only due of the 1-dimensional nature of the high field limit. The transverse corrections 
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however, have greatly changed. First we note an important additional term in \V±p\ 2 /p 
which does not depend on the strong field intensity. Perpendicular gradients are favored by 
this term. The next term in Vj_p, which is small compared to the constant term, is the same 
as found by Tomishima and Shinjo but with a factor +3. This term and the third term in 
Vj_p decrease with B Q . Therefore, contrary to Tomishima and Shinjo the total transverse 
corrections increase with Bq. This result modifies the transverse pinch. 

In summary, as expected, the behavior of the total corrections produces an anisotropic 
density profile. This anisotropy is caused by a reduction of dimensionality due to the high 
field limit. Our approach, which is technically simple, allows us to extend the calculus 
to finite temperatures and to include some Landau levels. But since the minimization of 
eqn. (^) is rather complicated, we leave its numerical solution for a forthcoming article. 
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